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Abstract: Sexual differentiation in Bombyx mori is controlled by sex-specific splicing of 
Bmdsx, which results in the omission of exons 3 and 4 in a male-specific manner. In 
B. mori, insulin-like growth factor II mRNA-binding protein (Imp) is a male-specific factor 
involved in male-specific splicing of Bmdsx. Male-specific Imp mRNA results from 
the male-specific inclusion of exon 8. To verify the link between histone methylation 
and alternative RNA processing in Imp, we examined the effects of RNAi-mediated 
knockdown of several histone methyltransferases on the sex-specific mRNA expression of 
Imp. As a result, male-specific expression of Imp mRNA was completely abolished when 
expression of the H3K79 methyltransferase DOT1L was repressed to <10% of that in 
control males. Chromatin immunoprecipitation-quantitative PCR analysis revealed a higher 
distribution of H3K79me2 in normal males than in normal females across Imp. RNA 
polymerase II (RNAP II) processivity assays indicated that RNAi knockdown of DOT1L in 
males caused a twofold decrease in RNAP II processivity compared to that in control 
males, with almost equivalent levels to those observed in normal females. Inhibition of 
RNAP Il-mediated elongation in male cells repressed the male-specific splicing of Imp. 
Our data suggest the possibility that H3K79me2 accumulation along Imp is associated 
with the male-specific alternative processing of Imp mRNA that results from increased 
RNAP II processivity. 
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1. Introduction 

Alternative splicing of pre-mRNA is an essential mechanism in the regulation of differential gene 
expression that can produce functionally distinct proteins from a single gene based on the 
developmental or physiological state of the cells in a multicellular organism [1,2]. Recent studies 
estimate that 90% of human genes are alternatively spliced [3], and several thousand different mRNA 
isoforms can be produced from a single gene. Although many examples describe how alternative 
splicing regulates gene expression, the mechanisms involved are less well understood [4-8]. 

Alternative splicing is thought to be regulated by the interaction of splicing factors and splicing 
enhancers (or silencers). Alternative splicing regulatory mechanisms have been investigated, structural 
models of spliceosomes have been proposed, and many RNA regulatory elements have been 
characterized; however, the emerging complexity of alternative splicing regulation suggests that these 
approaches have not sufficiently described how alternative splicing is regulated. Recent provocative 
studies point to a key function of chromatin structure and histone modification in alternative splicing 
regulation [9-11]. For example, the nucleosome occupancy level was lower in cassette exons than in 
constitutively spliced exons [12-14]. H3K36me3 and H3K9ac were related to the exon-skipping event 
of NCAM [15] and the levels of H3K36me3 differed in mutually exclusive exons of several genes 
among different cell types. Moreover, a genome-wide study across different species revealed that 
H3K36me3 was depleted in skipped exons [13,16]. Other histone modifications such as H3K4mel, 
H3K4me3, H3K27me3, and H3K9mel have been associated with the alternative splicing events of 
FGFR2 [17]. H3K4me3 was suggested to affect the alternative splicing events of CHD1 [18], while 
H3K9me3 was associated with multiple exon skipping of CD44 [19]. A recent genome- wide 
chromatin immunoprecipitation (ChlP)-seq analysis of histone H3 methylation in mammals revealed 
that alternative exons are preferentially marked with H3K4mel, H3K27me3, and H3K79me2, while 
being marked with H3K4me2, H3K4me3, and H3K36me3 at significantly lower levels [20]. 

In the silkworm Bombyx mori, the chromosomal sex determination mechanism is distinct from that 
of Drosophila melanogaster, with female (ZW) being the heterogametic sex and male (ZZ) the 
homogametic sex. The female sex in B. mori is determined by the presence of a dominant feminizing 
factor, Feminizer (Fern), on the W chromosome [21]. Note that no sex-specific regulatory Sxl homo log 
has been isolated from B. mori [22], and no tra homolog has been found in the Bombyx genome [23]. 
Despite these differences, a B. mori dsx homolog (Bmdsx) has been implicated in sex determination [24]. 
The primary transcript of the Bmdsx gene is alternatively spliced in males and females to yield 
sex-specific mRNAs that encode male-specific (BmDSXM) and female-specific (BmDSXF) 
polypeptides [25]. We found that unlike Drosophila dsx, the Bmdsx female exon is devoid of putative 
TRA/TRA-2 binding sites [25]. Instead, the splicing inhibitor BmPSI and a B. mori insulin-like growth 
factor II mRNA-binding protein (Imp) regulate male-specific splicing of Bmdsx [26,27]. Imp is 
localized on the Z chromosome and is expressed in a male-specific manner in various tissues. In male 
cells, the male-specific Imp mRNA is formed as a result of the inclusion of exon 8 and the promoter-distal 
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poly(A) site choice, whereas non-sex-specific polyadenylation occurs at the promoter-proximal 
poly(A) site downstream of exon 7 [28]. The molecular mechanisms underlying the sex-specific 
splicing regulation of this gene remain unclear. 

To verify the link between histone methylation and alternative RNA processing in Imp 
mRNA production, we investigated the effects of RNAi-mediated knockdown of several histone 
methyltransferases (HMTases) on sex-specific mRNA expression of Imp. Notably, the male-specific 
expression of Imp mRNA was completely abolished when expression of the H3K79 methyltransferase 
DOT1L was repressed to <10% of that in control males. Here, we provide several lines of evidence 
suggesting that H3K79me2 accumulation along Imp is associated with male-specific alternative RNA 
processing in Imp mRNA production, resulting from increased RNAP II processivity. To our 
knowledge, this is the first report to associate histone modification with the regulation of sex-specific 
alternative splicing. 

2. Results and Discussion 

2.1. Results 

2.1.1. Knockdown of DOT1L Abolished Male-Specific Expression of the Imp mRNA 

Recent genome-wide ChlP-seq analyses revealed that alternatively spliced exons are preferentially 
marked with H3K4mel, H3K27me3, and H3K79me2 [20]. Furthermore, a genome- wide study across 
different species revealed that H3K36me3 was depleted in skipped exons [13,16]. To investigate 
whether these epigenetic marks are associated with male-specific splicing of Imp pre-mRNA, we 
performed RNAi knockdown of several histone methyl transferases (HMTases) such as ASH2, EZH2, 
SETD2, and DOT1L known to modify H3K4, H3K27, H3K36, and H3K79, respectively, in embryos. 
Microinjection of dsRNA into B. mori embryos has been used successfully in many studies, although 
silencing levels vary [29]. siRNAs were injected into eggs during the early embryonic stage 6-8 h after 
oviposition, a developmental period known to be sensitive to RNAi-mediated gene knockdown [30]. 
Total RNA was extracted from each egg 4 days after injection. As shown in Figure 1A, qRT-PCR 
confirmed a significant reduction in EZH2, SETD2, and DOT1L transcript levels in embryos injected 
with siRNAs targeting these HMTase-coding genes. Injection of ASH2 siRNA failed to reduce the 
level of the target gene mRNA even though we used several siRNA sequences. Therefore, we focused 
on the knockdown effects of SETD2, DOT1L, and EZH2 on the expression of the male-specific Imp 
mRNA (Imp M ). RNAi knockdown of SETD2 and EZH2 had no influence on the expression of Imp M 
(Figure 1C, lanes 3, 4, 8, and 9). Notably, the expression of Imp M was completely abolished when the 
DOT1L expression level was repressed to <10% of that in control males (Figure 1C, lane 5). Five of 
six examined individuals whose DOT1L level was less than 10% also showed the disappearance of 
male-specific Imp expression. Further study is required to determine whether a similar effect on the 
expression of Imp M occurs when the expression levels of SETD2 and EZH2 are repressed <10% of that 
in control males. 
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Figure 1. The effect of histone methyltransferase (HMTase) knockdown on sex-specific 
splicing of Imp. (A) Quantification of HMTase mRNA expression 4 days after siRNA 
injections using quantitative reverse transcription (qRT)-PCR. Elongation factor 2 (EF-2) 
served as an internal standard. Error bar: SD; n = 8-24 individuals. ** p < 0.01, Student's 
?-test; (B) A schematic diagram of alternative splicing in Imp pre-mRNA. Exons are 
numbered and displayed as boxes. The gray box indicates the male-specific exon. The 
V-shaped lines above and below the diagram denote the splice variants observed in males 
and females. Imp contains two poly(A) sites. The proximal promoter site located within 
intron 7 is utilized in a non-sex-specific manner. The distal promoter site is selected in a 
male-specific manner and exists near the end of exon 8. The arrows indicate the 
approximate location of primers used for RT-PCR in C; (C) The male-specific Imp mRNA 
(Imp M ) was detected by RT-PCR and analyzed in a 1% agarose gel. The upper panel 
depicts expression of Imp M , and the lower panel shows the amplification of the Actin3 
transcript, which served as a positive control for RNA extraction and RT-PCR. Sex 
identification of each egg was performed by PCR amplification of the W-specific random 
amplified polymorphic DNA (RAPD) marker Rikishi. The expression levels of targeted 
HMTases relative to the negative control embryos in each individual examined are 
indicated below each lane. 
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2.1.2. DOT1L Knockdown Affects Male-Specific Splicing of Imp Pre-mRNA 



The above results indicate that DOT1L knockdown led to the loss of male-specific Imp expression 
in males. Two possible explanations may account for this phenomenon: DOT1L knockdown repressed 
Imp transcription or downregulation of DOT1L inhibited the splicing between exons 7 and 8 in Imp 
pre-mRNA. To examine these possibilities, we performed comparative analyses of Imp M and the Imp 
transcript common to both sexes {Imp ). RT-PCR analyses using primers described in Figure 2A 
demonstrated that DOT1L knockdown had little or no effect on the Imp mRNA expression (Figure 2B). 
In contrast, Imp M transcript was not detected in a male with the DOT1L expressed to <10% of that in a 
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control male (Figure 2B, lane 2). Next, using qRT-PCR, Imp mRNA expression was compared with 
that of the Imp c mRNA. As shown in Figure 2C, extremely low expression of Imp M mRNAs were 
observed in all males with DOT1L levels <10% of those in control males relative to Imp mRNA. 
These results indicate that DOT1L knockdown exclusively decreased the expression of the Imp M 
transcript. To rule out the possibility that down-regulation of DOT1L by RNAi affects Imp RNA 
transcription, qRT-PCR was performed to measure the Imp mRNA level that represents total Imp 
mRNA expression. As shown in Figure 2D, the Imp mRNA level in DOT1L knockdown males with 
DOT1L expression <10% of that in control males was decreased to a level less than one -third of that in 
control males (Figure 2D, left panel). In contrast, a more prominent reduction in Imp M mRNA was 
observed in the same DOT1L knockdown males (Figure 2D, right panel). These results indicate that 
DOT1L knockdown affects Imp RNA transcription but the male-specific splicing of Imp is repressed 
more severely by DOT1 depletion, leading to the loss of the Imp M expression. These results highlight 
the link between DOT1L and male-specific alternative splicing of Imp pre-mRNA. To investigate 
whether the decreased level of the Imp M mRNA caused by DOT1L depletion affects sex-specific 
splicing of Bmdsx, RT-PCR analysis was performed using primers designed to allow discrimination 
between female- and male-specific Bmdsx transcripts (Figure 2 A, lower panel). As shown in 
Figure 2E, decreased expression of Imp M induced the expression of the female-specific Bmdsx 
(Bmdsxfl and Bmdsxf2) mRNAs in addition to the male-specific Bmdsx (Bmdsx M ) mRNA (Figure 2E, 
lane 4). Only Bmdsxfl expression was observed males whose Imp M expression was severely 
diminished by DOT1L depletion (Figure 2E, lane 3). These results were consistent with our previous 
data showing that downregulation of Imp M by RNAi in male cells increases female-specific splicing of 
Bmdsx [27]. Sometimes a doublet band was seen in the Z)077Z-knockdown females (Figure 2E, lane 5). 
The upper band corresponded to Bmdsxfl and the lower band was expected to be a splice variant that 
lacks the third exon. At present we do not know exactly the reason for the appearance of the splice 
variant that lacks exon 3 in the DO 77Z-knockdown female. 

2.1.3. High Levels of H3K79me2 Favor Inclusion of Male-Specific Exon of Imp 

The above results support the possibility that H3K79 methylation marks are directly or indirectly 
associated with the regulation of the male-specific splicing of Imp pre-mRNA. We next performed 
comparative mapping of H3K79me2 across the alternatively spliced regions of Imp in females and 
males by ChlP-qPCR. As shown in Figure 3B, higher distributions of H3K79me2 were observed in 
males than in females across Imp (Figure 3B, left panel). Although the H3K79me2 enrichments were 
not limited to the alternatively spliced exon 8, most of the significant differences between males and 
females were observed in exon 8. In contrast, no significant differences in the levels of H3K79me2 
over Bmdsx were observed (Figure 3B, middle panel). To investigate whether the distribution of the 
other histone marks across Imp differed between males and females, we performed comparative 
mapping of H3K36me3 across Imp in females and males by ChlP-qPCR. As a result, no significant 
difference between males and females was observed in the levels of H3K36me3 over Imp (Figure 3B, 
right panel). To examine whether siRNA-mediated knockdown of DOT1L reduces the level of 
H3K79me2 accumulation, DOT1L siRNA- injected embryos were subjected to ChlP-qPCR analyses. 
Although higher distributions of H3K79me2 were observed in males than in females across Imp in 
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negative control individuals (Figure 3C, left panel), DOT1L knockdown exerted a marked influence on 
the accumulation of H3K79me2, with an at least 10-fold reduction throughout the regions examined in 
males (Figure 3C, right panel). Western blotting analysis demonstrated that DOT1L siRNA-injection 
efficiently reduced the total H3K79me2 level (Figure 4C). Since Dotl in yeast and its homologs 
appear to be solely responsible for H3K79 methylation [31], reduction in the accumulation of 
H3K79me2 across Imp by DOT1L depletion could be attributed to a reduction in the amount of 
DOT1L interacting with Imp. These results provide overwhelming support for the specificity of the 
ChIP data. However, reduced accumulation of H3K79me2 in females was less than that in males. This 
difference in the extent of reduction in H3K79me2 between males and females could be attributed to 
differences in £)077Z-knockdown efficiencies between males and females, as shown in Figure 1C and 
Figure 2E. Mean DOT1L expression relative to control individuals was 27.6% in DOT1L knockdown 
females and 15.1% in Z)077Z-knockdown males. At present we cannot explain the sexual difference in 
knockdown efficiencies. Although loss of H3K79me2 may affect indirectly the distribution of 
H3K4mel, H3K27me3, or H3K36me3, such indirect effect seems unlikely to affect the male-specific 
splicing of Imp since RNAi knockdown of ASH2, EZH2, and SETD2 had no influence on the 
expression pattern of Imp mRNA (Figure 1C). Overall we conclude that methylation of histone H3 
across Imp involves the HMTase DOT1L and high levels of H3K79me2 favor inclusion of exon 8. 

To verify the link between the high distributions of H3K79me2 and inclusion of alternative exons in 
genes other than Imp, we examined the effects of Z)077Z-knockdown on a Bombyx homolog of 
Sex-lethal gene {Bm-Sxl). The primary transcript of the Bm-Sxl gene is alternatively spliced to yield 
two isoforms, Bm-Sxl-L and Bm-Sxl-S [22]. Bm-Sxl-L consists of eight exons, while Bm-Sxl-S is a splice 
variant that lacks the second exon (Figure 4A). qRT-PCR analysis demonstrated that Bm-Sxl-L/ 'Bm-Sxl-S 
ratio was more than twofold higher in DO 77Z-knockdown embryos than in negative control embryos 
(Figure 4B). This result indicated that the £)077Z-knockdown relatively increased the inclusion of 
alternatively spliced exon (exon 2). In contrast, SisTZ^-knockdown caused no influence on the 
Bm-Sxl-L 'Bm-Sxl-S ratio. Western blotting analysis demonstrated that DOT1L siRNA- injection and 
SETD2 siRNA- injection efficiently reduced the total H3K79me2 level and total H3K36me3 level, 
respectively (Figure 4C). To investigate whether the high distributions of H3K79me2 is correlated 
with DOTlL-dependent inclusion of Bm-Sxl exon 2, we performed mapping of H3K79me2 around the 
exon 2 by ChlP-qPCR. The H3K79me2 mark was specifically enriched on the regions around exon 2 
as compared with the H3K36me3 mark (Figure 4D). These results suggest that enrichment of 
H3K79me2 could be correlated with DOTlL-dependent inclusion of alternative exons. 
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Figure 2. Effect of DOT1L knockdown on male-specific splicing of Imp. (A) The upper 
panel shows the location of the reverse transcription (RT)-PCR primers used for detection 
of Imp mRNA. The arrows above the diagram indicate the primers that were used for 
amplification of Imp M , and the arrows below the diagram show the primers used to amplify 
Imp mRNA transcribed from a region common to both sexes {Imp ). The lower panel 
indicates a schematic diagram of alternative splicing in the Bmdsx pre -mRNA. The gray 
boxes indicate the female-specific exons. The arrows point to the approximate locations of 
the primers used for RT-PCR in E; (B) Expression of Imp M and Imp c was detected by 
RT-PCR with primers illustrated in A and analyzed in a 1% agarose gel. The cDNA 
samples examined in lanes 1, 2, and 3 were identical to those used in lanes 1, 5, and 6, 
respectively, in Figure 1C. The bottom panel shows amplification of the ActinS transcript, 
which served as a positive control for RNA extraction and RT-PCR. The expression of 
DOT1L relative to the negative control embryos in each individual is indicated above each 
lane; (C) The ratio of Imp M to Imp c was analyzed by qRT-PCR; (D) Quantification of Imp 0 
(left panel) or Imp M mRNA expression (right panel) by qRT-PCR. The elongation factor 2 
(EF-2) served as an internal standard. SD; n = 5 individuals. * p < 0.05, *** p < 0.001, 
Student's Z-test; (E) Female- or male-specific splicing of Bmdsx pre-mRNA was detected 
by RT-PCR and analyzed on a 1% agarose gel. The upper panel shows expression of Imp M , 
and the middle panel indicates the female- and male-specific splicing products of Bmdsx 
(Bmdsxfl, Bmdsxf2 and Bmdsxm, respectively). The lower panel shows amplification of the 
Actin3 transcript, which served as a positive control for RNA extraction and RT-PCR. The 
expression levels of DOT1L and Imp M relative to those of the negative control embryos in 
each individual are indicated above each lane. 
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Figure 3. Increased H3K79 methylation across Imp in males. (A) Schematic representation 
of Imp and Bmdsx genes showing the distribution of quantitative (q) PCR amplicons used 
in the analysis; (B) Mapping of H3K79me2 across Imp (left panel) and Bmdsx (middle 
panel) and of H3K36me3 across Imp (right panel) in female (black) and male (white) larval 
tissues by chromatin immunoprecipitation (ChlP)-qPCR. Values represent the means ± SE 
of six qPCR values from one representative of five independent experiments. * p < 0.05, 
** p < 0.01, Student's West; (C) ChIP assays with antibodies to H3K79me2 and H3 and 
chromatin prepared from 60-pooled negative control embryos of each sex (left panel) or 
60-pooled DOT1L siRNA- injected embryos of each sex (right panel). The relative 
enrichment of H3K79me2 on EF-2 exon2 or along Imp was quantified by qPCR using 
primer sets indicated in A and expressed as a fraction of histone H3. Values represent the 
means ± SE of two independent qPCR assays from one representative of two independent 
experiments. * p < 0.05, Student's r-test. The percentage of input was normalized to 
unmodified H3. 
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Figure 4. Effect of DOT1L knockdown on alternative splicing of Bm-Sxl and distribution 
of H3K79me2 around the alternatively spliced exon in Bm-Sxl. (A) Schematic 
representation of Bm-Sxl showing the distribution of qPCR amplicons used in the analysis; 
(B) The ratio of Bm-Sxl-L to Bm-Sxl-S was analyzed by qRT-PCR. SD; n = 5 individuals. 
** p < 0.01, Student's £-test; (C) Western blotting analysis of H3K79me2 protein extracted 
from negative control or DOT1L siRNA- injected embryos using anti-H3K79me2 and 
anti-H3 antibodies (lower left panel). Quantification of H3K79me2 protein levels, as 
detected by Western blotting analysis (upper left panel). The intensity of each band was 
measured using Bioimage Analyser LAS 1000. H3K79me2 protein level was normalized to 
the H3 protein level. Values represent the means ± SE of six bands from one representative of 
two independent experiments. ** p < 0.01, Student's £-test. The same analysis was performed 
on H3K36me3 protein extracted from negative control or SETD2 siRNA-injected embryos 
using anti-H3K79me2 and anti-H3 antibodies (upper and lower right panels); (D) Mapping 
of H3K79me2 (black) and H3K36me3 (white) around Bm-Sxl exon 2 in larval tissues by 
ChlP-qPCR. Values represent the means ± SE of six qPCR values from one representative 
of five independent experiments. 

A 

Bm-Sxl-L 




2.1.4. Male and Female Differences in RNAP II Processivity in the Imp Gene 

Our data suggest a link between male-specific alternative splicing of Imp pre-mRNA and 
higher accumulation of H3K79me2. H3K79me2 modification is tightly associated with active 
transcription [32-34]. Moreover, several observations have suggested a close relationship between 
RNAP II dynamics and alternative splicing [35,36]. Based on these observations, we analyzed whether 
RNAP II processivity in Imp differs between males and females. To this end, we utilized an 
observation previously identified by others that transcription by a slower mutant RNAP II results in an 
increase in the ratio between promoter-proximal and promoter-distal pre-mRNA [15,36]. qRT-PCR 
of Imp pre-mRNA with primer sets located at each end of Imp (Figure 5A) showed that the 
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distal/proximal (D/P) pre-mRNA ratio was more than twofold higher in males than in females 
(Figure 5B). In contrast, no significant difference in the D/P ratio was observed in the control gene, 
B. mori elongation factor-1 (EF-1) (Figure 5C). These data suggest that a link between male-specific 
alternative splicing of Imp pre-mRNA and higher RNAP II processivity at this locus. As shown in 
Figure 5B, DOT1L knockdown reduced the D/P ratio in the male to a level similar to that in the normal 
female. Taken together with the ChlP-qPCR data indicated in Figure 3C, this result supports the 
possibility that higher distribution of H3K79me2 across Imp results in increased RNAP II processivity 
at this locus. 

2.1.5. Suppression of Male-Specific Imp Pre-mRNA Splicing by Inhibitors of Nucleotide Biosynthesis 

Sensitivity to MPA and the base analog 6AU are hallmarks of transcription elongation defects in 
yeast [37-39]. Both drugs cause the depletion of cellular nucleotide substrate pools required by RNA 
polymerases [40,41]. In vitro, RNAP II complexes pause and arrest more frequently under conditions 
of limiting amounts of nucleotides [42], suggesting that the in vivo hypersensitivity to these 
compounds is due to increased dependence on factors that promote elongation by RNAP II. Treatment 
with these drugs causes inhibition of exon skipping [36,43] because the extent of skipping of 
alternative exons correlates with the elongation rate of RNAP II [36]. Based on these reports, we 
investigated the effects of these inhibitors on male-specific splicing of Imp pre-mRNA. Male cultured 
cells (NIAS-Bm-Ml) that were not growth-arrested were treated with these inhibitors at concentrations 
established previously [43]. The efficiency of male-specific splicing was estimated by the ratio of 
Imp M /Imp c mRNA. The Imp M /Imp c ratio in male cells treated with each drug relative to that in the 
negative control cells is indicated in Figure 5D. In the MPA experiment, an approximately threefold 
decrease in the Imp M IImp c ratio was observed. A similarly high reduction in the Imp M IImp c ratio was 
observed when the male cells were treated with 6AU. Importantly, no significant difference was 
detected in the level of total Imp mRNA between control cells and inhibitor-treated cells, indicating 
that the reduction in the Imp M /Imp c ratio represented simply the reduction in male-specific Imp mRNA 
expression. Taken together with measurements of RNAP II processivity in Figure 5B, the simplest 
explanation for these results is that drug-induced reduction in transcript elongation led to the inhibition 
of the male-specific splicing of Imp pre-mRNA. 

2.1.6. Embryonic Lethality Caused by DOT1L Knockdown 

Above results indicated that DOT1L knockdown caused decreased expression of Imp M , leading to 
the expression of the female-specific Bmdsx mRNA in male embryos (Figure 2E). To evaluate whether 
DOT1L does indeed play an important role in sexual differentiation, we investigated the effect of 
DOT1L knockdown on the development of sexual phenotypes. The highest hatch rate of the control 
siRNA- injected embryos in six trials was 21.7% (Table 1), which was still lower than that reported by 
another group [30]. Presumably, this difference was caused by technical issues related to microinjection. 
Similarly low hatchability in negative control dsRNA- injected eggs is reported by the other group [44]. 
Compared with the control embryos, nearly all the embryos injected with DOT1L siRNA did not hatch, 
suggesting embryonic lethality (Table 1). One hatched larva was obtained from a male egg injected 
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with DOT1L siRNA that survived to the adult stage. This male had normal fertility and its genital 
organs showed no abnormalities when viewed under a dissecting microscope (data not shown). 

Figure 5. Male and female differences in RNAP II processivity in Imp. (A) Schematic 
diagram of Imp and EF-1 showing the distribution of the proximal and distal amplicons 
(black bars) used for qPCR analysis; (B,C) RNAP II processivity was determined as a ratio 
of the proximal and distal pre-mRNA sequences (distal/proximal) of Imp (B) and EF-1 (C). 
The abundance of each pre-mRNA was determined by quantitative (q)PCR. Values 
represent the means ± SE of six qPCR values from three individuals; (D) Inhibitors of 
nucleotide biosynthesis suppress male-specific splicing of Imp pre-mRNA. qRT-PCR 
analysis was performed to calculate the ratio of Imp M to Imp c . The Imp M IImp ratio in male 
cells treated with 0.1 mM mycophenolic acid (MP A) or with 0.2 mM 6-azauracil (6AU) 
(C) is relative to that in the negative control cells in each experiment. Values represent the 
means ± SE from three individual experiments. * p < 0.05, ** p < 0.01, Student's £-test. 
Values below each graph indicate total Imp mRNA expression in the cells examined 
relative to those in the negative control cells. 
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Table 1. Effects of DOTIL siRNA injection on egg development. 



siRNA 


Injected 

eggs 


Early or mid-stage 
embryonic lethal 


Late-stage 
embryonic lethal 


Viable 


N. C. siRNA 


115 


58 (50.4%) 


32 (27.8%) 


25 (21.7%) 


DOTIL siRNA 


263 


188(71.5%) 


74 (28.1%) 


1 (0.4%) (male) 



2.2. Discussion 



Imp in B. mori has been identified as a male-specific RNA-binding protein involved in the 
regulation of male-specific splicing of Bmdsx [27]. The pre-mRNA of Imp undergoes sex-specific 
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RNA processing. In male cells, male-specific Imp mRNA contains exon 8 and the distal promoter 
poly(A) site choice, whereas non-sex-specific polyadenylation occurs at the proximal promoter 
poly(A) site downstream of exon 7. Here, we found that DOT1L knockdown affects male-specific 
splicing of Imp pre-mRNA, leading to loss of male-specific Imp expression (Figure 2B,C). In support 
of this result, higher distributions of H3K79me2 were observed in males than in females or in DOT1L 
knockdown males across Imp (Figure 3B,C). Strong link between enrichment of H3K79me2 and 
DOTlL-dependent inclusion of alternative exons was also observed in Bm-Sxl (Figure 4). Comparative 
analysis of RNAP II processivity indicated that higher distribution of H3K79me2 across Imp was 
correlated to increased RNAP II processivity at this locus (Figure 5B). Inhibition experiments using 
inhibitors of RNAP II elongation suggested that the higher elongation rate was tightly associated with 
male-specific RNA processing of Imp pre -mRNA (Figure 5D). Together, our data suggest that 
H3K79me2 accumulation along Imp is associated with male-specific alternative RNA processing in 
Imp mRNA production, resulting from increased RNAP II processivity. 

The alternative RNA processing pattern of Imp pre -mRNA closely resembles that found in 
Drosophila polo pre-mRNA. Polo, which is a cell cycle gene, also contains a proximal and a distal 
poly (A) site in the 3 'untranslated region (UTR) to produce alternative mRNA that differ in their 3'UTR 
length [45]. In a mutant Drosophila strain that displays a reduced RNAP II elongation rate, RNAP II 
occupancy along polo is altered and the proximal poly(A) site is used 3.5-fold more efficiently than in 
wild-type flies [46]. An increase in proximal poly(A) site usage was also observed in five other 
alternatively polyadenylated transcripts in Drosophila. These results in Drosophila show that the 
kinetics of RNAP II can determine alternative poly(A) site selection. As shown in Figure 5D, 
significant reduction in the Imp M IImp c ratio was observed when the male cells were treated with 
inhibitors of RNAP II elongation. This result indicates that drug-induced repression of transcript 
elongation leads to the reduction in utilization of distal poly(A) site, resulting in the relative increase in 
proximal poly(A) site usage. Slow RNAP II presumably exposes the proximal poly(A) site on the 
nascent transcript to the polyadenylation machinery for a longer time before RNAP II transcribes the 
distal poly(A) site [46]. Therefore, the proximal poly(A) site is processed before the distal poly(A) site 
is transcribed, suggests a "first come, first served" mechanism. This resembles the extra domain I 
(EDI) alternative splicing mechanism described previously, whereby a slow RNAP II preferentially 
included the normally excluded alternative EDI exon because it allowed ample assembly time for the 
spliceosome machinery [36,47]. 

Eleven-nineteen lysine -rich leukemia gene (ELL) family proteins are essential components of the 
super elongation complex (SEC) and increase the catalytic rate of transcription elongation by RNA 
polymerase II [48-50]. ELL2 knockdown by siRNA affects the alternative pre-mRNA processing of 
the immunoglobulin heavy chain (IgH) gene, which is accompanied by reduced H3K4 and H3K79 
methylation [50]. In addition to the ELLs, the SEC contains the MLL translocation partners 
AF4/FMR2 family member 1 (AFF1; also known as AF4), AFF4, eleven-nineteen leukemia (ENL) 
and ALL 1 -fused gene from chromosome 9 (AF9) [47]. In this complex, ENL is linked, not only with 
all members of the AF4 protein family that occur as MLL fusion partners, but also with pTEFb and 
DOT1L [51,52]. Notably, several frequent MLL fusion partners seem to coordinate DOT1L activity 
with a protein complex that stimulates the elongation phase of transcription by phosphorylating the 
carboxy-terminal repeat domain of RNA polymerase II [49,50]. 
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Figure 6. Model for the regulation of alternative splicing of Imp pre-mRNA by H3K79me2 
and DOT1L. (A) In males, a transcription elongation factor (X) promotes the formation of 
SEC on Imp, causing H3K79 methylation by DOT1L. Higher RNAP II processivity due to 
H3K79me2 accumulation does not allow enough time to complete 3'-end processing at the 
proximal poly(A) site (AAGAAA), leading to exclusive use of the distal 3' splice site of 
exon 8 and the distal poly(A) site (AAUAAA). While the distal poly(A) signal sequence is 
perfectly matched to the most canonical poly(A) signal hexamer, the proximal poly(A) 
signal sequence is consistent with a single-nucleotide variant of the canonical hexamer. 
Therefore, the polyadenylation machinery prefers the distal poly(A) site rather than the 
proximal poly(A) site when both poly(A) sites are present on the nascent transcript 
simultaneously. The protein product from the male-specific Imp transcript induces the 
male-specific splicing of Bmdsx pre-mRNA; (B) In females, the presence of a dominant 
feminizing factor, Fern, on the W chromosome directly or indirectly represses the 
expression of X, leading to failure of SEC formation. Decreased accumulation of 
H3K79me2 caused by the loss of SEC slows RNAP II processivity, providing sufficient 
time to recruit cleavage factors for 3'-end processing at the proximal poly(A) site. The absence 
of the male-specific Imp induces the female-specific splicing of Bmdsx pre-mRNA. 



Male-saecific splicing 




fast 



B 

Female-specific splicing 




H3K79me2 modification is tightly associated with active transcription [32,33]. Milcarek et al. 
speculated that conversion of monomethylated H3K79 into di- and trimethylated forms is correlated 
with the transition from low to high level gene transcription, due most likely to a decrease in the 
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histone-DNA interaction [53]. Because H3K79 resides within the histone core, its methylation may 
facilitate DNA unwinding from the histone, allowing the downstream chromatin to open more readily 
and be transcribed more efficiently. 

Based on these previous findings, we propose a possible model for the regulatory mechanism 
underlying sex-specific alternative splicing of Imp pre-mRNA (Figure 6). In male cells, exclusive 
expression of a transcription elongation factor (X), such as ELL family proteins, promotes the 
formation of SEC on Imp, causing H3K79 methylation by DOT1L. Higher RNAP II processivity due 
to H3K79me2 accumulation does not allow enough time to complete 3'-end processing at the 
proximal, non-sex-specific poly(A) site, leading to exclusive use of the 3' splice site of exon 8. The 
male-specific Imp protein induces the male-specific splicing of Bmdsx pre-mRNA (Figure 6A). We 
have found recently that male-specific Imp bound immediately downstream of the proximal poly(A) 
site and promoted male-specific splicing of its pre-mRNA [28]. Therefore, after the male-specific Imp 
protein has been produced, the protein product may inhibit use of the proximal poly(A) site and promote 
the splicing of intron 7, leading to exclusive use of the 3' splice site of exon 8. In female cells, Fern 
directly or indirectly represses the expression of X, leading to failure of SEC formation. Decreased 
accumulation of H3K79me2 caused by the loss of SEC slows RNAP II processivity, providing 
sufficient time to recruit cleavage factors such as CF I, CF II, and/or poly(A) polymerase (PAP) for 
3'-end processing at the proximal poly(A) site. The absence of the male-specific Imp induces the 
female-specific splicing of Bmdsx pre-mRNA (Figure 6B). 

In this model, both the proximal and distal poly(A) sites are exposed to the polyadenylation factors 
at the same time in male cells. How then is the distal poly(A) site exclusively selected? While the 
distal poly(A) signal sequence is perfectly matched to the most canonical poly(A) signal hexamer 
AAUAAA, the proximal poly(A) signal sequence is AAGAAA, which is consistent with a single-nucleotide 
variant of the canonical hexamer. The motifs that are functional in vertebrates are AAUAAA and its 
highly conserved variants (e.g., AUUAAA, UAUAAA, AGUAAA, AAGAAA) [54]. Among those 
hexamers, the canonical AAUAAA was reported to be present in 53% of the mRNAs; in contrast, a 
single-nucleotide variant AAGAAA was found in only 3% of the mRNA [55]. Therefore, one can 
reasonably suppose that the polyadenylation machinery prefers the distal poly(A) site rather than the 
proximal poly(A) site when both poly(A) sites are present on the nascent transcript simultaneously. 
This scenario is consistent with the kinetic coupling model for the regulation of alternative splicing by 
RNAP II elongation [56]. In this model, when a proximal suboptimal (weak) 3' splice site and a 
downstream canonical (strong) 3' splice site are presented simultaneously to the splicing machinery, 
the strong 3' splice site could easily outcompete the weak site, resulting in alternative exon skipping. 
Our results indicate that both alternative polyadenylation and alternative splicing depend on RNAP II 
kinetics. In the present study, we were unable to assess precisely whether DOT1L plays a crucial role 
in regulating sex determination or sexual differentiation of B. mori because nearly all the embryos 
injected with DOT1L siRNA died before hatching. Consistent with our results, Z)oZ7Z-deficient 
embryos died between 9.5-10.5 days post coitum due to developmental abnormalities, including 
growth impairment, angiogenesis defects in the yolk sac, and cardiac dilation [57]. In D. melanogaster, 
grappa {gpp) is an ortholog of DotlL [58]. gpp is an essential gene identified in a genetic screen for 
dominant suppressors of pairing-dependent silencing where a Poly comb-group (Pc-G)-mediated 
silencing mechanism necessary for the maintenance of parasegment identity during embryo 
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development [59]. As is the case in D. melanogaster, DOT1L in B. mori may be required to maintain 
developmental gene expression through Pc-G-mediated mechanism. 

Assuming that our model presented above is valid, an ELL family protein — such as ELL2 — might 
be a key factor in facilitating male differentiation as a result of inducing the male-specific splicing of 
Imp pre-mRNA. Recently, we found that an ELL2 homolog is expressed in embryos in early 
developmental stages. Ongoing investigations are aimed at determining whether the ELL2 homolog is 
involved in regulating the sex-specific splicing of Imp pre-mRNAs. 

3. Experimental Section 

3.1. Silkworm Strains 

The Bombyx mori non-diapausal and white egg strain (pnd-wl) was kindly provided by Kenichi 
Moto of RIKEN (Wako, Osaka, Japan). The S-2 strain, in which the females have the T (W; 2, 5) p B + re 
(black egg, black larvae) genotype and the males have the + pB , re (red egg, white larvae) genotype, 
was established in our laboratory. The former strain was used primarily for RNAi experiments and the 
latter for ChlP-quantitative PCR (qPCR) and gene expression analysis. The developing eggs were 
enclosed in a plastic case and incubated at 25 ± 2 °C with sufficient humidity. Larvae were reared on 
an artificial diet (Nihon Nosan, Yokohama, Japan) at 25 ± 2 °C. 

3.2. Preparation of siRNAs 

cDNA sequences predicted to encode ASH2, EZH2, SETD2, and DOT1L were retrieved by 
tBLASTn searches of the KAIKObase (http://sgp.dna.affrc.go.jp/KAIKObase/) using human ASH2, 
EZH2, SETD2, and D0T1L as query sequences (Figure 7). siRNA targeted to four HMTases (ASH2, 
EZH2, SETD2, D0T1L) were designed as described previously [60]. All sequences used in RNAi 
experiments are listed in Table 2. Each siRNA was synthesized using the custom select siRNA 
synthesis service provided by Ambion (Austin, TX, USA). Silencer negative control #1 siRNA 
(Ambion, Austin, TX, USA) was used as a negative control in the siRNA experiments. 

3.3. Injection of siRNAs into Eggs 

siRNAs were injected into eggs as described previously [60]. Negative control siRNA, EZH2 
siRNA, and ASH2 siRNA were injected at 50 uM, and SETD2 siRNA and DOT1L siRNA were 
injected at 50 or 100 uM. 

3.4. Extraction of Total RNA and Genomic DNA 

Total RNA was extracted from each egg and from NIAS-Bm-Ml cells using Isogen (Nippon Gene, 
Tokyo, Japan) according to the protocol described previously [60]. Genomic DNA was recovered by 
ethanol precipitation from the intermediate and organic phases obtained in the RNA extraction process. 
The precipitated DNA was purified using the SimplePrep" reagent for DNA (Takara, Kyoto, Japan) 
according to the manufacturer's instructions. To perform molecular sexing of each egg, PCR was 
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performed using primers specific to the W chromosome genomic sequence according to the protocol 
described previously [60]. 

Figure 7. Predicted amino acid sequences of Bombyx mori ASH2, EZH2, SETD2, and 
DOT1L. (A) Alignment of Bombyx mori ASH2 (KAIKObase China Gene Model Gene No. 
BGIBMGA008025), Drosophila melanogaster ASH2, and human ASH2. The putative 
SPRY domain is denoted by horizontal line above the amino acid alignment; (B) Alignment 
of Bombyx mori EZH (KAIKObase EST clone ID: FS904534); (C) Drosophila melanogaster 
EZH, and human EZH; (C) Alignment of Bombyx mori SETD2 (KAIKObase China Gene 
Model Gene No. BGIBMGA003106), Drosophila melanogaster SETD2, and human 
SETD2. The putative AWS and SET domains are denoted by horizontal lines above the 
amino acid alignment; (D) Alignment of Bombyx mori DOT1L (KAIKObase FLcDNA 
clone name: ffbm34A09), Drosophila melanogaster Gpp (DOT1L ortholog), and human 
DOT1L. The putative DOT1 domain is denoted by horizontal line above the amino acid 
alignment. Amino acid identity is denoted by black boxes. 
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Table 2. siRNA sequences used in RNAi experiments. 



Target gene siRNA 



Sense 



Antisense 



SETD2 



DOT1L 



EZH2 



ASH2 



393 UGCCAGCUCUGAGUCUGAUUCAAU 
490 CAGUGUAGCUCAAGAGAUATT 

243 UUCCAAAGCAACUACAGAAUCGAUG 
353 AUUUACUCGCUUUACUUUGTT 

38 GACAACCCAACAGGUACCAAUAAGA 

224 CGACGGGAAAGUGCAUGGUGAUAAA 

157 GACCGGCCUCUAGUCAAAUUCAAGA 

167 UAGUCAAAUUCAAGAGCCACCUGUA 



AUUGAAUCAGACUCAGAGCUGGCA 

UAUCUCUUGAGCUACACUGTT 
CAUCGAUUCUGUAGUUGCUUUGGAA 

CAAAGUAAAGCGAGUAAAUTT 
UCUUAUUGGUACCUGUUGGGUUGUC 
UUUAUCACCAUGCACUUUCCCGUCG 
UCUUGAAUUUGACUAGAGGCCGUC 
UACAGGUGGCUCUUGAAUUUGACUA 



3.5. Reverse Transcription (RT)-PCR Analyses 

RT-PCR reactions performed according to the protocol described previously [60]. The primer 
sequences and PCR conditions utilized in this study are indicated in Supplementary Table 3. 



Table 3. Primer sequences and PCR conditions utilized in this study. 



Target gene 


Primers 


Sequence (5' to 3') 


Denature 


Annealing 


Elongation 


N° cycles 




FDSX-F2 


CGCCTTACCGCAGACAGGCAG 


98 °C 


57 °C 


57 °C 




Bmdsx 












35 




FDSX-R4 


GCGCAGTGTCGTCGCTACAAGG 


10 s 


30 s 


60s 






BmlMPFl 


ATGGACGGTGACATGTCTCAAG 


98 °C 


55 °C 


57 °C 




Imp M 












35 


BmlMPRl 


TCATCCCGCCTCAGACGATTG 


10 s 


30 s 


90s 






IMPE4F1 


TCCCATAATAATCTCATTGGAC 


98 °C 


55 °C 


57 °C 




Imp c 












35 


IMPE7R1 


AATGTGAACGGTGGTCTCGTG 


10 s 


30 s 


90s 






BA3F1 


AGATGACCCAGATCATGTTCG 


98 °C 


57 °C 


57 °C 




Actin3 












26 




BASR1 


GAGATCCACATCTGTTGGAAG 


10 s 


30 s 


30 s 






BmSxlFl 


ATTAATCATCATAAAGCTACG 


98 °C 


57 °C 


57 °C 




Bm-Sxl 












35 




BmSxlRl 


AATCCGTAACTGTAGCCAGTC 


10 s 


30 s 


30 s 





3.6. Quantitative Real-Time RT-PCR (qRT-PCR) 

qRT-PCR assays were performed according to the protocol described previously [60]. All primer 
sequences used in this study are listed in Table 4. The BmEF-2Fl and BmEF-2Rl primers were used 
to amplify elongation factor-2 (EF-2) as an internal standard for quantification [61]. 



Table 4. Sequences of primers used for qRT-PCR. 



Target gene 


Primers 


Sequence (5' to 3') 


SETD2 


SETD2qPCRFl 


CCTACAGGACATCTGGAGTTAC 


SETD2qPCRRl 


GAATCAGTACCAGCATTTAGATG 


DOT1L 


dotlqPCRFl 


AGAATCCGAACGACTCGACAG 


dotlqPCRRl 


CTGTTCTTGGTCTTCGTTCAAC 


EZH2 


EZH2F1 


GGTGTAGTGACAACCCAACAG 


EZH2R1 


TCTTAACTCCTGAGCTGTTCC 


ASH2 


ASHF2 


GGGGACCAGGTTCCACGAGTC 


ASHR2 


TACAGGTGGCTCTTGAATTTG 
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Table 4. Cont. 



Target gene 


Primers 


Sequence (5' to 3') 


Bmdsx 


BmdsxproFl 


TGCATGTTTCTTATTAATCAGCTAG 


promoter 


BmdsxproRl 


rn AAA rnrnrn x - ^ m A A A A S~\ s~~*\ : i i x - "i * x - ^ x - ^ * 

GTAAATTTCGTAAAAGCTGACCAG 


Bmdsx El 


ChlPExlF 


CCTGTACCACCAGTGGTGAAG 


ChlPExlR 


CTGACGGCGGTGGAGCGTATG 


Bmdsx 12 


x"~1 1 i t— v w r^t T " 

CWPI2F 


CTACATGAACAGTACCAGTCAG 


ChIPI2R 


rn A A S~\ rn A x - 1 AAA x - ^ m AAA rn A i x - ^ s~-i mm x - ^ 

GTAAGTACAAACTAAATAGCGTTC 


Bmdsx E3 


ChIPEx3F 


GTCGACGAGTACGCGAGGAAG 


ChIPEx3R 


m s~\ m i a m x - ^ a m s~\ rr~i a m/~lm x - ^ m/~l x - ^ x - ^ 

TGTGATGCATGTATCTGTCGC 


Bmdsx 13 


CWPI3F 


m A A m s~\ A x - ^ /imm x - mm x - ^ *^irxn A A m 

GTAACTGACCTTCTTGCTAATC 


CMPI3R 


s~-\ m s~\ m x - ^ x - ^ a / i v ■ v ■ v ■ i * mm A A m A m x - ^ x - ^ m x - ^ 

CTGTGCCATTTTATTAATATCGTC 


Bmdsx E4 


ChIPEx4F 


A m A m A A /^(m/^ y~<rny~(m A /~im/^rn/nrnrn/~i 

ATATAAGTGGTGTACTGTCTTC 


ChlPEx4R 


CCATAGATCCAATGTTACGAC 


Bmdsx 14 


CWPI4F 


GTTCAAACACATCGAAGCTAC 


ChIPI4R 


x - ^ m x - ^ x - ^ x - ^ a i a m A x~-i a x - ^ m i x - ^ x - ^ y^imm 

GTCCGAGATAGACTGGCCTTG 


Bmdsx E6 


ChIPEx6F 


x - ^ x~-i x - ^ A X - ^ A 1 X - ^ X - ^ y—A y—A s~\ A X - ^ A A S~\ m A A 

GGCACAGCGCCGACAAGTAAG 


ChIPEx6R 


A mrn/~<rn/^im/~<rn a s~\ a m a mrnin/^rn/~( a m 

ATTGTCTGTAGATATTCGTGATC 


BmIMP 


IMrprorz 


T"~T A A A TTT A A T'T' ATA A A A z - ^ A -p/^ 

1 1 AACjLA 111 AA 1 1 A 1 AACjAACjA 1 L 


promoter IS 


IMPproR2 


CTAGAATCTGCGATTACATAC 


BmIMP VIS 


IMPE1F1 


m s~ ^ x - ^ mm x - ^ a s~\ ' i ' A x - ^ m x - ^ x - ^ /-im A m A x - ^ 

TCCGTTCAGTACTCGCTATAC 


IMPE1R1 


m s~ i mm A x - ^ x - ^ m A m s~ ^ m x - ^ a m A i a mm x - ^ 

TCTTACCTATCGTCATAGATTC 


BmIMP US 


IMPI1F2 


A mmm i m A A A A m A m x - ^ m s~ ^ m A m x - i 

ATTTGGTAAAATAGTCTCGTATC 


IMPHR2 


A x - ^ mm ^ m i a m A s~ ^ i x - ^ x - ^ mm A A 

ACCTTGTGATACGGGGTTAAC 


BmIMP prolL 


IMPElLproFl 


x - ^ x~- ^ m x - ^ A x - ^ y~ ^ x - ^ mmm AAA x - ^ x - ^ x - ^ 

GCTGCCCCACCCTTTAAACCG 


IMPElLproRl 


CTCGATCGTGCTGACTCTAGC 


BmIMP EIL 


IMPEIFI 


TTTCAAGTATACTCCTTCTATAG 


IMPE1R1 


mm x - ^ x - ^ x - ^ s~~i a mmmm x - ^ * x~-* x~ ^ a x~-i a mm 

TTCGCCATTTTGAGCAGATTG 


BmIMP III, 


IMPI1F3 


x - ^ AAA m x - ^ a s~ ^ A m A mm x - ^ mm x - ^ 

CAAATGGGCACATATTGTTGG 


IMPI1R3 


/"(rnmrn A A /"( >nrnrnm/nynrn/~( A m/~t 

GTTTAAGCGCTTTCGTGATGG 


BmIMP E7 


IMPE7F1 


ATGCGGGAAGAAGGTTTTATG 


IMPE7R1 


AATGTGAACGGTGGTCTCGTG 


flm/MP 17 


IMPI7F1 


GTGCATAAATCCACAGAACAG 


IMPI7R1 


mm A /^m/^l A s~"\ AAA x - ^ m x - ^ a x~-i a a i m A x - ^ 

TTACTCAGAAACTCAGAAGTAC 


Sw/MP E8 


IMPE8F1 


CGTCTGAGGCGGGATGAGAAC 


IMPE8R1 


m AAA mm x - ^ x - ^ x - ^ x - ^ x - ^ x - ^ a a m x - ^ a s~\ s~-\ a x~-i 

TAAATTCGCCGCAATCAGCAG 


£7*7 promoter 


EF-lproFl 
EF-lproRl 


TATATCAATTTTGGTGCAAGAATGG 

m A A m A A m A mm ✓—i m A mmi^m A m x - ^ x - ^ a x - i x - ^ x - ^ 

GTAATAATATTCTATTCTATCCACCG 


TiF-7 E2 


EF-1E2F1 
EF-1E2R1 


TGGCGATGGAGGCGGAGAAG 

x~~< rp ^-"1 » » /~ir-r-tr-r-t/—i A /^/^rp/-^m < nm/-^^( 

CTCAACTTCCCAGCTGTCTGC 


£F-7 13 


EF-1I3F1 
EF-1I3R1 


ACTTACTTATTTATGATCATGCGTC 
GCTAACCACAATTATATTTGTGGAG 


E6 


EF-1E6F1 
EF-1E6R1 


TACAGGTCATTTCTGCACGTAAG 
TCATCCCAGTTAACTGTTGGATC 


EF-1 111 


EF-1I11F1 
EF-1I11R1 


TTCATGGACTACATTTTACCTTGG 
CTAAGCTCTTCTAAAAGAGATGAGC 
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Table 4. Cont. 



Target gene 


Primers 


Sequence (5 to 3 ) 


EF-1 112 


EF-1I12F1 
EF-1H2R1 


GCATTAATATTAATTCCACCACAAG 
CACACCTCACTGCTCTTCCGC 


Bm-Sxl 


SxlilFl 


GGCTAAACTATCTTCAACAAG 


promoter 


SxlilRl 


CGGTCACCGTTCTCGTGAAAG 


Bm-Sxl El 


SxlelFl 
SxlelRl 


GCCAGTCCAAATGGACGAATC 

s~\ mm y - ^ a y^imi y-^ A i^lmmm y - ^ A s~\ m y~"^ A s~~^ 

GTTCACTGACTTTCGAGTGAG 


Sm-Sil E2 


Sxle2Fl 
Sxle2Rl 


GCCTACTCGAACAATAAAAAAG 

mm s~-\ AAA y-^ A A mm y-^ AAA y - ^ m y~ *i y^imi y-^ 

TTCCAAAGAATTGAAACTCCTG 


Bm-Sxl 12 


Sxli2Fl 
Sxli2Rl 


m/^s A A m/~l A y~s A A a m/^lm/^1 A mmm/^l y-^ 

TGAATCAGAACATCTCATTTGG 

s~-\ y—A a a y~s y - ^ x - ^ /-Nm y -- * y - ^ /-im a y - ^ x~~"i m a a y - ^ 

CCAAGCCGCTGCCTACCTAAC 


Sm-Sx/ E3 


Sxle3Fl 
Sxle3Rl 


CGAGGCAGAGCGGGTTCGAAC 

x~""i /~imm y - ^ a m y~ *i A /~im/^l y-^ A y - ^ a y-^ /~im/^l m y - ^ 

CCTTCATCACTCGACAGCTCTC 


BmIMP 


IMPE2F1 


ATCCTCAAAGGTACTCATCAG 


Proximal 


IMPE2I2R 


y~^ y - ^ a mi y~s y - ^ a m y - ^ a /~im/^l A A y - ^ a a m A 

GCATGCATCACTCAACAATAC 


BmIMP Distal 


IMPE7F2 


GGATCATCGGCAAAGGCGGAC 


IMPI7R2 


A y* - ^ y - ^ a /"1mm y-^ y-^ A m y~ *i A mm y - ^ a m A y - ^ y—A 

AGCACTTGGATCATTCATACC 


£"F-7 Proximal 


EF-1E2F1 
EF-1E2I2R 


TGGCGATGGAGGCGGAGAAGG 

y* - ^ A A A A A A A S~-\ AAA S~\ X - ^ A mm X - ^ A T~l y - "* X - ^ 

GAAAAAGAAGAAAGCATTCATGC 


EF-1 Distal 


EF-1I11E12F 
EF-1E12R2 


A m y-^ A m A /im m A mm A A x - ^ : i i x~ ^ a mm x - ^ 

ATGATACTGTATTAACTGCATTC 

mm x - ^ a y~s x - ^ a mi mi mi mi y~s A y - "* A x - ^ r ■ 1 x~-* 

TTCAGGATTTTGAGACCCTGG 


Sm/MP E7-E8 


IMPE7F1 
BmlMPRl 


A m y-^ x - ^ y - ^ x - ^ y - ^ a A y~S A A y-^ y - ^ mmmm A m y-^ 

ATGCGGGAAGAAGGTTTTATG 

mi x - ^ a mi y~ *i x - ^ y - ^ y - ^ y - ^ y - ^ m x - ^ a s~\ a x - ^ y - ^ a mm 

TCATCCCGCCTCAGACGATTG 


5W/MPE1-E3 


BmlMPFl 
IMPE3R1 


ATGGACGGTGACATGTCTCAAG 

x~~"i a m x - ^ x - ^ a mm x - ^ a a x - ^ x - ^ x - ^ x^ mmm A m y~^ 

CATCCATTCAACCCGTTTATG 


BmIMP VIS 


BmlMPFl 


ATGGACGGTGACATGTCTCAAG 


IMPE2R2 


y~* /~i /^(m/~( inrp/~1rp/^ y~l A /— i mi /— i mi /— i y— ( A A y— ( 

GCCTGCTCTGGACTCTCGAAG 


BmIMP E1L 


IMPE1LF1 


AATCTGCTCAAAATGGCGAAG 


IMPE2R2 


GCCTGCTCTGGACTCTCGAAG 


Sw-&c/ S 


BmSxlF2 
BmSxlSRl 


ACTCGCGTTACCTATTTAAC 
GTACTGCTGTTGGATTTGGTC 


Bm-Sxl L 


BmSxlF2 
BmSxlLRl 


ACTCGCGTTACCTATTTAAC 
CTGCTGTTGGATTTGATTTTC 



3.7. ChIP Experiments 

The fifth instar day-3 larvae fat bodies were treated with 1 x cold phosphate-buffered saline (PBS) 
containing 1% protease inhibitor cocktail (Roche, Basel, Switzerland). To cross-link samples, a 
formaldehyde solution (50 mM HEPES-KOH [pH 7.5], 100 mM NaCl, 1 mM EDTA [pH 8.0], 0.5 mM 
EGTA [pH 8.0], 11% formaldehyde) was added to a final concentration of 1% for 15 min at room 
temperature. Eggs were mashed with a pestle in the formaldehyde solution to prepare samples for ChIP 
experiments. The reaction was stopped by adding a glycine solution to a final concentration of 
125 mM for 5 min at room temperature and washed twice with lx cold PBS. The fixed fat bodies were 
disrupted with a Polytron (KINEMATIC A, Lucern, Switzerland) (10,000 rpm, 20 s) and centrifuged at 
2800 rpm for 5 min at 4 °C. The samples were resuspended in lysis buffer (10 mM Tris-HCl [pH 8.0], 



Int. J. Mol. Sci. 2014, 15 



6791 



100 mM NaCl, 1 mM EDTA [pH 8.0], 0.5 mM EGTA [pH 8.0], 0.1% Na-deoxycholate, 
0.5% iV-lauroylsarcosine, 1% protease inhibitor) and sonicated (pulsed 30 s, paused 1 min x 7 sets) on 
ice to avoid overheating. Triton-X was added to the resulting lysate at a final concentration of 1% and 
centrifuged at 15,000 rpm for 10 min at 4 °C. An aliquot (200 jxL) of each sonicated sample served as 
an input DNA control. For ChIP, 300 uL of the supernatant was incubated with Dynabeads Protein G 
(Invitrogen, Carlsbad, CA, USA) with the appropriate antibody overnight at 4 °C. The beads were 
washed twice with lx cold PBS. Next, 5 ug of antibodies were added in 300-uL blocking solution 
(0.5% bovine serum albumin [BSA]/PBS) and incubated overnight at 4 °C. The antibodies used in this 
study were anti-H3 antibody (abl791, Abeam, Cambridge, UK), anti-H3K36me3 antibody (ab9050, 
Abeam, Cambridge, UK), anti-H3K79me2 antibody (ab3594, Abeam, Cambridge, UK), and anti-rabbit 
IgG antibody (12-370, Millipore, Billerica, MA, USA). The beads were washed once with low salt 
buffer (20 mM Tris-HCl [pH 8.0], 150 mM NaCl, 2 mM EDTA [pH 8.0], 0.1% SDS, 1% Triton X-100), 
twice with high salt buffer (20 mM Tris-HCl [pH 8.0], 400 mM NaCl, 2 mM EDTA [pH 8.0], 0.1% SDS, 
1% Triton X-100), five times with RIPA buffer (50 mM HEPES-KOH [pH 7.5], 500 mM LiCl, 1 mM 
EDTA [pH 8.0], 1% NP-40, 0.7% Na-deoxycholate), and once with TE (Tris + EDTA) + 50 mM 
NaCl. To elute the histone-DNA complex, elution buffer (50 mM Tris-HCl [pH 8.0], 10 mM EDTA 
[pH 8.0], 1% SDS) was added to the beads and incubated at 65 °C for 15 min followed by 
centrifugation to obtain the supernatant. The supernatants were further incubated at 65 °C for 6 h, and 
then incubated with RNase A (Roche, Basel, Switzerland) at 37 °C for 2 h followed by incubation with 
proteinase K (20 mg/mL, Takara, Kyoto, Japan) at 55 °C for 2 h. Finally, the DNA was extracted with 
phenol chloroform and diluted in ddH 2 0. To calculate the amount of target sequence in the 
immunoprecipitated chromatin, we performed real-time qPCR as described above. 

3.8. Mycophenolic Acid (MP A) and 6-Azauracil (6AU) Treatments 

NIAS-Bm-Ml cells were maintained in IPL-41 (Invitrogen, Carlsbad, CA, USA) with 10% fetal 
bovine serum (FBS) (HyClone, Thermo Scientific, Waltham, MA, USA) under a humidifying 
atmosphere at 26 ± 1 °C and 2 mL of cultured cells were seeded onto a 35-mm dish. For MPA 
treatment, 2.5 uL of MPA (16 mg/mL in dimethyl sulfoxide (DMSO); Wako, Osaka, Japan) were 
added to each dish. For the 6AU treatment, 2.5 uL of 6AU (2 mg/mL in DMSO; Wako, Osaka, Japan) 
were added to each dish. The concentration of each inhibitor was determined as reported 
previously [43]. Total RNA was isolated at 1 or 3 days after treatment, respectively, according to the 
protocol described above. 

4. Conclusions 

Male-specific alternative splicing of Imp pre-mRNA was repressed by DOT1L depletion in male 
embryos. Consistent with this finding, higher distributions of H3K79me2 were observed in males than 
in females across Imp. Comparative analysis of RNAP II processivity indicated that RNAP II 
processivity was higher in males than in females at this locus. Inhibition experiments using inhibitors 
of RNAP II elongation suggested that the higher elongation rate was closely associated with 
male-specific splicing of Imp pre-mRNA. Taken together, our data suggest that greater accumulation 
of H3K79me2 along Imp in males causes increased RNAP II processivity, leading to male-specific 
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alternative RNA processing in Imp mRNA production. Furthermore, knockdown of DOT1L caused 
embryonic lethality. 
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